Revised version March 6, 2000 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



SUPERNOVAE VERSUS NEUTRON STAR MERGERS AS THE MAJOR EPROCESS SOURCES 

Y.-Z. Qian 

School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455; qian@physics.umn.edu 

Revised version March 6, 2000 



o 
o 
o 



> ■ 

<n : 

(N ■ 

m ■ 

o : 
o . 
o ■ 

^ ' 

i ' 
o ■ 

c3 



ABSTRACT 

I show that recent observations of r-process abundances in metal-poor stars are difficult to explain 
if neutron star mergers (NSMs) are the major r-process sources. In contrast, such observations and 
meteoritic data on 182 Hf and 129 I in the early solar system support a self-consistent picture of r-process 
enrichment by supernovae (SNe). While further theoretical studies of r-process production and enrich- 
ment are needed for both SNe and NSMs, I emphasize two possible direct observational tests of the SN 
r-process model: gamma rays from decay of r-process nuclei in SN remnants and surface contamination 
of the companion by SN r-process ejecta in binaries. 

Subject headings: Galaxy: evolution — nuclear reactions, nucleosynthesis, abundances — supernovae: 
general — stars: neutron 



1. INTRODUCTION 

All of the actinides and approximately half of the stable 
nuclei with mass numbers A > 100 in the solar system 
were produced by rapid neutron capture, the r-process. 
While the astrophysical site of this process remains a 
mystery, the main candidate environments are neutrino- 
heated ejecta from core-collapse supernovae (SNe; Woosley 
& Hoffman 1992; Meyer et al. 1992; Takahashi, Witti, & 
Janka 1994; Woosley ct al. 1994) and decompressed ejecta 
from neutron star mergers (NSMs; Lattimer & Schramm 
1974, 1976; Symbalisty & Schramm 1982; Freiburghaus, 
Rosswog, & Thielemann 1999b). To be a major source 
for the r-process, an environment must satisfy two crite- 
rions: one on reproducing the solar r-process abundance 
pattern and the other on supplying the total amount of r- 
process material in the present Galaxy. Using more phys- 
ical parametrizations than the previous approach based 
on constant neutron number density and temperature 
(e.g., Kratz et al. 1993), several recent studies (Hoffman, 
Woosley, & Qian 1997; Meyer & Brown 1997; Freiburghaus 
et al. 1999a) derived the r-process conditions that can pro- 
duce relative abundance patterns with peaks at A ~ 130 
and 195 as observed in the solar system. At present, it 
seems that conditions in models of neutrino-heated ejecta 
from SNe are deficient (e.g., Qian & Woosley 1996) while 
those in models of decompressed ejecta from NSMs are 
promising (e.g., Freiburghaus et al. 1999b) for an re- 
process. However, due to the uncertainties in the theoreti- 
cal models (e.g., Qian & Woosley 1996; Freiburghaus et al. 
1999b), a reliable comparison of the actual conditions in 
these two environments with the r-process conditions can- 
not be made yet to help establish or discriminate either 
environment as the r-process site. Furthermore, for both 
environments the models showed that if the ejecta were 
composed of r-process material, then the amount provided 
by a single event combined with the estimated number of 
SNe or NSMs over Galactic history would be adequate to 
account for the present Galactic r-process inventory (e.g., 
Mathews & Cowan 1990; Qian & Woosley 1996; Rosswog 
et al. 1999). Therefore, the above two criterions can- 
not readily identify SNe or NSMs as the major r-process 



sources. 

In this Letter I discuss a phcnomenological approach to 
test SNe and NSMs as the major r-process sources. By 
considering mixing of the ejecta from an individual event 
with the interstellar medium (ISM), I show that observa- 
tions of metal-poor stars are difficult to explain if NSMs 
are the major r-process sources (§2). I further show that 
a self-consistent picture of r-process enrichment by SNe 
based on the same consideration is supported by mete- 
oritic data on 182 Hf and 129 I in the early solar system and 
by observations of metal-poor stars (§3). To emphasize 
the importance of observations in establishing the major 
r-process sources, I conclude with a discussion of two pos- 
sible direct tests of the SN r-process model: gamma rays 
from decay of r-process nuclei in SN remnants and surface 
contamination of the companion by SN r-process ejecta in 
binaries (§4). 

2. ^-PROCESS ABUNDANCES AND NEUTRON STAR 
MERGERS 

I assume that as far as the r-process is concerned, the 
mergering of a neutron star (NS) with a black hole (BH) 
is similar to that of two neutron stars. These two kinds 
of events are collectively referred to as NSMs. The rate 
of such events in the Galaxy is quite uncertain. Here I 
take an average rate (/q SM ) ~ (10 4 yr) 1 over Galac- 
tic history, which is at the very high end of various esti- 
mates (e.g., Phinney 1991; Bethe & Brown 1998; Arzou- 
manian, Cordes, & Wasserman 1999). Numerical simula- 
tions of a NS-NS merger event by Rosswog et al. (1999) 
showed that a total mass M^ SM ~4x (10- 3 -10- 2 )M© 
of decompressed NS material may be ejected. Then the 
grand total from all the past NSMs over Galactic history 
of t G ~ 10 10 yr is M| SM </£ SM }t G _ 4 x (i()3-4O 4 )M . 
This is roughly equal to the present Galactic r-process in- 
ventory X%*M G ~ 10 4 M Q , where X^ r ~ 10" 7 is the 
total r-process mass fraction of nuclei with A ;> 100 in 
the solar system (Kappeler, Beer, & Wisshak 1989) and 
Mq ~ 10 11 Mq is the total Galactic mass in stars and 
gas. So it seems that NSMs could be the major r-process 
sources. 
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However, the amount of ejecta from a single NSM event 
discussed above is too much to explain the observed r- 
process abundances in metal-poor stars. This can be seen 
by considering mixing of the ejecta from each event with 
the ISM. Rosswog et al. (1999) showed that the total en- 
ergy of the NSM ejecta is at most comparable to the SN 
explosion energy (~ 10 51 erg). Therefore, when its en- 
ergy/momentum is dispersed in the ISM, this ejecta can 
mix with at most the same amount of material as swept 
up by a SN remnant, M mix «3x 1O 4 M (e.g., Thorn- 
ton et al. 1998). Consequently, if all of this ejecta were 
r-process material as required to account for the present 
Galactic r-process inventory, an ISM enriched by a single 
NSM event would have a total r-process mass fraction 



^C,NSM — 1-3 x 10 7 



M NSM 



4 X lO" 3 Mr. 



'3 x 1O 4 M 



(1) 

The solar r-process mass fractions of elements with 100 <J 

100<A<130 



A < 130 and A > 130 are X 
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^0,r 



4 x 10~ 8 (Kappeler et al. 1989). According to equa- 
tion (|j) , X*jij SM would be approximately equal to X@f r ~ 

100<A<130 .. „„,. 

X Q r + X^> 130 even for the lowest M^ SM of inter- 
est. Therefore, whether the NSM ejecta were composed 
of r-process elements with 100 <J A < 130 or A > 130, 
or both groups in solar proportion, equation (Q) predicts 
that abundance ratios of e.g., Ag (A ~ 107) and/or Eu 
(A ~ 153) with respect to hydrogen in an ISM enriched 
by a single event would be at least comparable to the corre- 
sponding solar r-process values (Ag/H) Q)T . and (Eu/H) Q)T .. 
This predicted level of r-process enrichment is in disagree- 
ment with recent observations of r-process abundances in 
metal-poor stars, as Ag/H in stars with [Fe/H] ~ —1.3 to 
—2.2 are ~ 10-10 2 times lower than (Ag/H)©.,, (Crawford 

et al. 1998) while Eu/H in stars with [Fe/H] 3 are 

- 30-10 3 lower than (Eu/H) , r (McWilliam et al. 1995; 
Sneden et al. 1996, 1998). 

The following interpretation of equation (|Tj) gives some 
insights into how NSMs would explain the present Galac- 
tic r-process inventory and may help appreciate why this 
explanation is disfavored by observations. The Galaxy 
can be divided into ~ 3 x 10 6 units each having a mass 
Mmjx ~ 3 x 10 4 A/q. This division has a physical mean- 
ing as M m i x is the maximum mass within which the ejecta 
from an individual NSM event can be distributed. Due 
to the rather low rate of NSMs in the Galaxy, on average 
at most one NSM event occurred in a unit over Galac- 
tic history. Therefore, in order to account for the present 
Galactic r-process inventory, each unit would have to be 
enriched with an approximately solar r-process mass frac- 
tion by a single NSM event. As shown above, this picture 
of r-process enrichment is inconsistent with the observed 
r-process abundances in metal-poor stars. Furthermore, 
due to the high rate of SNe in the Galaxy, many SNe oc- 
curred in a unit where a single NSM event also occurred 
at sometime in Galactic history (see §4). As Fe enrich- 
ment of this unit was provided by these SNe, stars formed 
at different times in this unit would have varying [Fe/H] 
but either zero or approximately solar r-process mass frac- 
tion if NSMs were the major r-process sources. This is in 
disagreement with the observed correlation between abun- 



dances of r-process elements and Fe at [Fe/H] J> —2.5 
(Gratton & Sneden 1994; Crawford et al. 1998; see also 
McWilliam et al. 1995). 

Of course, the above discussion of mixing of the NSM 
ejecta and the ISM is oversimplified. For example, one 
could imagine that a smaller than average fraction of r- 
process ejecta was mixed into the ISM near the edge of a 
NSM remnant. In this case stars formed near the edge of 
the remnant would have r-process mass fractions smaller 
than that given by equation (|IJ). However, one would also 
expect that less than average enrichment was not unduly 
pervasive and a significant fraction of the metal-poor stars 
would have the r-process enrichment indicated by equation 
(0). The fact that no such stars have been observed sug- 
gests a difficulty of the NSM r-process model in explaining 
the observations at low metallicities. Furthermore, even 
if the observed r-process abundances in metal-poor stars 
could be attributed to pervasive less than average enrich- 
ment by NSMs, one would still face the other difficulty in 
explaining the observed correlation between abundances 
of r-process elements and Fe at [Fe/H] ^ — 2.5. As Fe 
enrichment was controlled by SNe occurring at a much 
higher rate, widely- varying degrees of mixing of the r- 
process ejecta in an already existing NSM remnant with 
Fe produced by fresh SNe would result in large scatter in 
r-process abundances over a broad range of [Fe/H]. This is 
in disagreement with the rather early onset of the correla- 
tion between abundances of r-process elements and Fe at 
[Fe/H] rj — 2.5. In summary, observations of metal-poor 
stars would be difficult to explain if NSMs were the major 
r-process sources. 

3. .R-PROCESS ENRICHMENT BY SUPERNOVAE 

In contrast to the case of NSMs, observations of metal- 
poor stars as well as meteoritic data on 182 Hf and 129 I in 
the early solar system support a self-consistent picture of 
r-process enrichment by SNe (Qian & Wasserburg 2000, 
hereafter QW; Wasserburg & Qian 2000; hereafter WQ). 
In this picture the ejecta from each SN event is mixed with 
an average mass M m ; x w 3 x 10 4 M Q of ISM swept up by 
the SN remnant (e.g., Thornton et al. 1998). It is as- 
sumed that the SN rate per unit mass of gas /Q N /M gas is 
approximately constant over Galactic history (this seems 
reasonable as the star formation rate is proportional to the 
gas mass). Consequently, an average ISM in the Galaxy 
is enriched by newly-synthesized material from SNe at a 
frequency 



SN 



fSN 
JG 

x M 



(10 7 yr)- 



3 x 10 4 Mi 
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fSN 
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(30 yr)" 



10 10 Aid 

Mtras 



, (2) 



where /,| N /Mg as is estimated using quantities for the 
present Galaxy. 

Meteoritic data on 182 Hf (with a lifetime fi82 = 1-3 x 
10 7 yr) and 129 I (f 129 = 2.3 x 10 7 yr) in the early so- 
lar system shed important light on the r-process and 
its association with SNe. Wasserburg, Busso, & Gallino 
(1996; see also QW) showed that the abundance ratio 
( 182 Hf/ 180 Hf) SS F = 2.4 x 10" 4 (Harper & Jacobsen 1996; 
Lee & Halliday 1995, 1998) is consistent with common 
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SNe injecting Hf into the ISM at a high frequency 
f H « ~ (10 7 yr) -1 over a uniform production time 
Tup ~ 10 yr before solar system formation (SSF). How- 
ever, the abundance ratio ( 129 I/ 127 I)ssf = 10~ 4 (Reynolds 
1960; see also Brazzle et al. 1999) must be accounted for 
by a different type of SNe occurring at a low frequency 
fc ~ /w/10 ~ (10 s yr)" 1 (Wasserburg et al. 1996; QW). 
Therefore, the meteoritic data require at least two dis- 
tinct types of SN r-process events: the high-frequency TL 
events producing heavy nuclei with A > 130 including 
182 Hf and the low-frequency C events producing light nu- 
clei with A < 130 including 129 1. The r-process produc- 
tion in the SN environments associated with the TL and 
C events was discussed in some detail by Qian, Vogel, & 
Wasserburg (1998a). 

The above picture of r-process production and enrich- 
ment by SNe has clear predictions for r-process abun- 
dances resulting from a single event (QW; WQ). For ex- 
ample, with an average ISM enriched by the TL events at 
a frequency fn ~ (10 7 yr) -1 , the solar r-process abun- 
dances of elements with A > 130 such as Eu were pro- 
vided by /-hTup ~ 10 3 TL events. This requires that the 
Eu abundance in an ISM enriched by a single TL event 
must be 

loge„(Eu) = loge , r (Eu) - log(/ w T UP ) ~ -2.5, (3) 

where the loge notation is defined as loge(Eu) = 
log(Eu/H) + 12 with Eu/H being the abundance ratio of 
Eu to hydrogen, and loge , r (Eu) = 0.51 is the value for 
solar r-process Eu (Kappeler et al. 1989). Similarly, the 
solar r-process abundances of elements with A < 130 such 
as Ag were provided by fcT\jp ~ 10 2 C events. This re- 
quires that the Ag abundance in an ISM enriched by a 
single C event must be 

loge £ (Ag) = log e0 , r (Ag) - log(/ £ T UP ) ~ -0.8, (4) 

where loge , r (Ag) = 1.19 (Kappeler et al. 1989) is used. 
The predictions in equations da) and (^) are in good agree- 
ment with observations of very metal-poor stars which 
were formed when only a small number of SNe had con- 
tributed r-process elements to the ISM. The observed 
loge(Eu) values for stars with [Fe/H] ~ —3 range from 
-2.5 to -0.9 (McWilliam et al. 1995; Sneden et al. 1996, 
1998), which can be accounted for by ~ 1-40 TL events 
with logew(Eu) ~ -2.5 from a single event (QW; WQ). 
In addition, Ag abundances at the level indicated by equa- 
tion (|) were observed in HD 2665 and BD +37° 1458 with 

[Fe/H] 2 (Crawford et al. 1998) and in CS 22892-052 

with [Fe/H] 3 (Cowan et al. 1999). 

In the above picture of r-process enrichment by SNe the 
total mass of r-process elements ejected in an TL event must 
be A^> 130 M mix /(/ w Tup) - lO~ 6 Af , while that in an 

100<A<130 , , 

C event must be ^ , r ~ M mix /(/ c Tup) - 1O- 5 M . 
A total ~ 1O~ 6 -1O -5 M of material can be naturally 
provided by the neutrino-heated ejecta from the proto- 
neutron star in a SN over a period ~ 1-10 s (e.g., Qian 
& Woosley 1996). However, whether this neutrino- heated 
ejecta is composed of r-process material is yet to be shown. 
The difference by a factor ~ 10 in the total amount of 
r-process ejecta between the TL and C events has been 
speculated to indicate that neutrino emission and hence, 



ejection of r-process material are terminated by the tran- 
sition of the proto-neutron star into a BH in the TL events 
while prolonged ejection is sustained by neutrino emission 
from a stable NS in the C events (Qian et al. 1998a). 

In summary, despite the lack of a complete model for 
successful r-process production in SNe, there is a self- 
consistent picture of r-process enrichment by SNe that can 
account for the meteoritic data on 182 Hf and 129 I in the 
early solar system and the observed r-process abundances 
in metal-poor stars. Furthermore, the observed correla- 
tion between abundances of r-process elements and Fe at 
[Fe/H] £ -2.5 (Gratton & Sneden 1994; Crawford et al. 
1998; see also McWilliam et al. 1995) can be understood as 
the result of sufficient mixing of r-process products and Fe 
from multiple SN events in this picture. In fact, the same 
picture was used by Wasserburg & Qian (WQ) as the ba- 
sic framework to explain the dispersion in abundances of 
heavy r-process elements such as Eu at [Fe/H] ~ —3. 

4. DISCUSSION AND CONCLUSIONS 

The amount of ISM to mix with the ejecta from an in- 
dividual NSM event is limited by the total energy of the 
event. On the other hand, due to the very low rate of 
NSMs in the Galaxy, a large amount of r-process ejecta 
would be required from each event to account for the 
present Galactic r-process inventory. When mixed with 
the ISM, this required amount of ejecta would result in 
abundances of r-process elements with A < 130 (such as 
Ag) and A > 130 (such as Eu) that are much too high (by 
factors ~ 10-10 2 for Ag and ~ 30-10 3 for Eu) compared 
with the observed values in metal-poor stars. Further- 
more, an average ISM received the ejecta from only ~ 1 
NSM event over Galactic history. If r-process enrichment 
of the ISM was provided by NSMs in this way while Fe 
enrichment was provided by many SNe, there would be 
no correlation between abundances of r-process elements 
and Fe, in disagreement with the observed correlation at 
[Fe/H] ^> —2.5. While the complexities in mixing of the 
ejecta with the ISM can affect the above considerations in 
detail, it is unlikely that they can remove the difficulty of 
the NSM r-process model in explaining the observations 
of metal-poor stars (especially the rather early onset of 
the correlation between abundances of r-process elements 
and Fe at [Fe/H] ps —2.5). Nevertheless, future numerical 
studies of r-process enrichment by NSMs accompanied by 
Fe enrichment by SNe should be interesting to pursue and 
may give a more definitive answer. 

In contrast, a self-consistent picture of r-process enrich- 
ment by SNe can be obtained by considering mixing of the 
ejecta from an individual event with the ISM. Here an av- 
erage ISM is enriched in r-process elements with A > 130 
by the TL events at a frequency fu ~ (10 7 yr) -1 and 
in those with A < 130 by the C events at a frequency 
fc ~ (10 8 yr) -1 . This picture can account for the mete- 
oritic data on 182 Hf and 129 I in the early solar system and 
the observed r-process abundances in metal-poor stars. 
Furthermore, sufficient mixing of r-process products and 
Fe from multiple SN events in this picture would result in 
the observed correlation between abundances of r-process 
elements and Fe at [Fe/H] ;> —2.5. The same picture was 
also used by Wasserburg & Qian (WQ) as the basic frame- 
work to explain the dispersion in abundances of heavy r- 
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process elements such as Eu at [Fe/H] ~ —3. However, a 
complete model of r-process and Fe production in the TL 
and C events is still lacking and should be investigated by 
future theoretical studies. 

In discussing r-process enrichment by NSMs I have as- 
sumed that the maximum amount of ISM to mix with the 
ejecta from an individual event is the same as swept up 
by a SN remnant. This is because the total energy of the 
NSM ejecta seen in numerical simulations (Rosswog et al. 
1999) is at most comparable to the SN explosion energy 
(~ 10 51 erg). For given conditions of the ISM, the expan- 
sion/evolution of the ejecta is essentially governed by its 
total energy. The large difference in the amount of ejecta 
between a NSM and a SN has no significant effect here as in 
both cases the mass of the swept-up ISM soon overwhelms 
that of the ejecta while the total energy remains more or 
less the same. I note that a small amount (<; 10~ 5 M Q ) 
of material might be ejected in highly-relativistic jets in a 
NS-BH merger event (Janka et al. 1999). However, the 
total energy of these jets is <J 10 51 erg (Janka et al. 1999) 
and their existence would not increase significantly the 
amount of ISM that could mix with the entire ejecta from 
the event. It takes ~ 10 6 yr for the energy (~ 10 51 erg) and 
the associated momentum of the ejecta to be dispersed in 
the ISM (e.g., Thornton et al. 1998), where the next NSM 
or SN would occur on a much longer timescale (~ 10 10 yr 
for NSMs and ~ 10 7 yr for SNe). This leaves substan- 
tial time for mixing of the ejecta with the swept-up ISM. 
However, the details of the mixing process require further 
studies. 

If, as argued here, SNe are the major sources for the 
r-process, then there are two possible direct observational 
tests: gamma rays from decay of r-process nuclei in SN 
remnants and surface contamination of the companion by 
SN r-process ejecta in binaries. Qian et al. (1998b, 1999) 
discussed gamma-ray signals from decay of long-lived re- 
process nuclei (with lifetimes ~ 10 3 -10 5 yr) in a nearby 
SN remnant and from decay of short-lived r-process nu- 
clei (with lifetimes ~ 1-10 yr) produced in a Galactic 
SN that may occur in the future. The nuclide 126 Sn is 
of particular interest (Qian et al. 1998b) as its lifetime 
(~ 10 5 yr) is much longer than the age (~ 10 4 yr) of the 
Vela SN remnant at a distance w 250 pc. In the picture 
of r-process enrichment by SNe discussed here (see also 
QW; WQ) the solar r-process mass fraction of nuclei with 
A < 130 resulted from ~ 10 2 C events (see §3). So a total 
mass Xg= 126 M mix /10 2 - 5 x 10- y M Q of 126 Sn nuclei are 
produced in each £ event, where Xq= 126 w 1.6 x 10~ 9 



(Kappeler et al. 1989) is the solar r-process mass frac- 
tion of 126 Tc, the stable daughter of 126 Sn. If the SN 
associated with the Vela remnant was an C event, then 
decay of 126 Sn in this remnant would produce gamma- 
ray fluxes ~ 10~ 7 7 cm~ 2 s _1 at energies E 7 = 415, 666, 
and 695 keV. Detection of these fluxes would require fu- 
ture gamma-ray experiments such as the proposed Ad- 
vanced Telescope for High Energy Nuclear Astrophysics 
(ATHENA, Kurfess 1994). As the Vela remnant contains 
a pulsar, such detection would also provide evidence for 
the speculated association between C events and SNe pro- 
ducing neutron stars (Qian et al. 1998a). 

The other test mentioned above takes advantage of the 
occurrence of SNe in binaries. The r-process ejecta from 
the SN would contaminate the surface of its binary com- 
panion. Some binaries would survive the SN explosion 
and acquire a NS or a BH in place of the SN progenitor. 
Therefore, an ordinary star observed to be the binary com- 
panion of a NS or a BH might show r-process abundance 
anomalies on the surface. To estimate the plausible level 
of such anomalies, I assume that a fraction ~ 10 ~ 3 of the 
entire SN ejecta (~ 10 Mq, mostly non-r-process material) 
would be intercepted by a low mass (~ 1 Mq) companion 
and then mixed with ~ 10~ 2 Mq of the surface material. 
If the SN was an 7i event, ~ 10~ 9 M Q of r-process ele- 
ments with A > 130 would be intercepted, while for an 
C event ~ 10~ 8 M Q of r-process elements with A < 130 
would be intercepted (see §3). These quantities arc to 
be compared with ~4x 10~ 10 Mq of the corresponding re- 
process elements in ~ 10 -2 Mq of the surface material in a 
companion star of solar r-process composition. So a SN in 
a binary could enhance significantly the surface r-process 
abundances in the companion star, especially if the SN 
was an C event. In view of the large overabundance of O, 
Mg, Si, and S recently observed in the companion star of a 
probable BH (Israelian et al. 1999), detection of r-process 
enhancement in similar binary systems seems promising. 
Such detection may also test directly the speculation by 
Qian et al. (1998a) that H events are associated with SNe 
producing BHs while C events are associated with SNe 
producing neutron stars. 
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